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Abstract 



\Q ' I review recent progress in simplicial quantum gravity in three and four 

dimensions, in particular new results on the phase structure of modified models 
of dynamical triangulations, the application of a strong-coupling expansion, 
and the benefits provided by including degenerate triangulations. In addition, 
I describe some recent numerical and analytical results on anisotropic crystalline 
£f) , membranes. 



O . 1 Simplicial Quantum Gravity in D > 2 

The task of formulating a consistent theory of quantum gravity in four dimensions, 
hence unifying the two pillars of modern physics, general relativity and quantum 
mechanics, is a formidable one and still unresolved. Many different approaches have 
been tried (see e.g. Ref. [[!]]), some have led to a deeper understanding of the problem 
but, as of yet, none can claim much success. In this article I will review the status, 
and recent progress, of one such approach, namely the attempt to make sense of an 
Euclidean path-integral quantization of general relativity via a discretization known 
as dynamical triangulations. 

In Euclidean quantum gravity a path-integral is written as a formal sum over 
D-dimensional Euclidean geometries (metrics g^ v ) 

V ^ e -s E [g^} m 

Vol(diff) { } 

weighted by the Einstein-Hilbert action 

i r JD , 



SbW = J d WM(-i* + 2A). (2) 

G is the Newton's constant, A the cosmological constant and R the scalar curvature. 
However, this path-integral formulation faces some sever problems and unresolved 
questions: 

• the theory is non-renormalizablc, 

• the action is unbounded from below, 
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• adding higher-order curvature terms leads to renormalizable, albeit non- 
unitary, theory, 

• how can Lorentzian signature be recovered from the Euclidean formulation. 

A proposed interpretation is to view Eq. (|2|) as the first part of an effective 
action, valid in the infrared limit, and that fine-tuning of the parameters in the 
full effective theory might reveal a non-trivial ultraviolet stable fixed point where 
a consistent quantum theory of gravity might be defined Q (a behavior akin to 
that of the non- linear c-model). While it is true that non-trivial fixed points in 
four-dimensional field theories are rare, some interesting results obtained in 2+e 
dimensions actually lend support to this idea Q. 

To address this proposal, a non-perturbative investigation of the theory is needed. 
One such is provided by a discretization of the path-integral Eq. (|l|) , usually follow- 
ing a prescription due to Regge . He proposed to approximate a smooth manifold 
with a given metric by a piecewise linear manifold (a triangulation) , on which both 
parallel transport and the integral of curvature are well-defined. Choosing a suffi- 
ciently fine triangulation should provide a good approximation to any continuum 
manifold. 

This idea has been applied to lattice gravity following one of two complimentary 
prescriptions, differing in how they implement the dynamical nature of the space- 
time metric: 

(a) Quantum Regge calculus 

Use a piecewise linear manifold with fixed connectivity but with varying edge 
lengths II . 

(b) Dynamical Triangulations 

Summing over triangulations with varying connectivity but fixed edge lengths || . 

In this review I will focus on the latter approach, for a recent review of quantum 
Regge calculus see e.g. Ref. M. 

1.1 Dynamical triangulations 

In models of dynamical triangulations the continuum integral over diffeomorphism 
inequivalent metrics is replaced by a discrete sum over all possible decompositions 
(gluings) of equilateral D-simplexes along their (D—l) -dimensional faces into sim- 
plicial manifolds, while requering that the neighbourhood of each vertex is a D-ball 
(excluding pseudo-manifolds) [^, |j| . Fixing all edge lengths to a constant a intro- 
duces a "reparametization invariant" cut-off in the model. 

In this approach the Einstein-Hilbert action is approximated by a simple function 
of the numbers iV, of i-dimensional simplexes on a given triangulations T. However, 
as these numbers are related through the Dehn-Sommerville and Eulers relations, the 
action can be taken to depend on only two of those numbers. The (grand-canonical) 
partition function becomes: 

Z^,k)=J2^ Nd E clT) eKND - 2 > W 
N D TeT Nn { > 
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Figure 1: The (p , g)-moves in three dimensions: (a) inserting/deleting a vertex, (6) replac- 
ing a triangle by an edge and vice verse. 

where [i and n represent the discrete cosmological and inverse Newtons's constants, 
and C(T) is the order of the automorphishm group of T — the number of equivalent 
labelings of the vertexes. 

The second sum is over a suitable ensemble T of D-dimensional triangulations 
(of volume Nd). A priori, different ensembles can be used provided they lead to a 
well-defined partition function. While this leaves ample freedom in defining different 
ensembles, the topology of the triangulations must be fixed else the partition func- 
tion is divergent in Np. So far, most simulations in D > 2 have used triangulations 
of spherical topology. 

The main justification for the dynamical triangulation approach comes from 
two-dimensions where variety of exact solutions exist, complemented by extensive 
numerical results, and which agree with the solutions of continuum Liouville theory.^ 

1.2 Simulation methods 

The partition function Eq. (J3j) is evaluated numerically using a Monte Carlo algo- 
rithm. The space of all triangulations is explored using a set of D+l local (topology 
preserving) geometric moves — the (p , q)— moves [jll|. In a (p,q)~move, where 
p = D + 1 — q, a (q — l)-subsimplex is replaced by its "dual" (p — l)-subsimplex, 
provided no manifold constraint is violated. The moves are ergodic in D < 4, i.e. 
any triangulation can be transformed into any other by an appropriate sequence of 
moves. An example of the (p,q)-moves in three dimensions is shown in n 1. 

In D > 2 fluctuations in the volume are necessary for an ergodic updating 
procedure and it is customary to simulate the quasi-canonical partition function: 

1 In contrast, simulations of quantum Regge calculus in two dimensions fail to reproduce the 
known continuum results [ [j"o[ . 
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= E e ~" ND £ ^e*-^-^, (4) 



where the volume fluctuates around a target volume iVo and the quadratic potential 
term added to the action ensures that, for an appropriate choice of 5, the fluctuations 
are small. 

To explore the phase structure of model Eq. (Q) powerful methods for probing 
the nature of the quantum geometry have been developed. The fractal structure of 
the triangulations is labeled by a set of critical exponents, e.g: 

(i) The string susceptibility exponent 7^ defines the singular behavior of the grand- 
canonical partition function: Z(\i) ~ Z rcg + (fi — /i c ) 2_7s . This in turn implies that 
the canonical partition function behaves asymptotically as 

Z Nd (k) ~ e^ N ° iV£ (K) - 3 , Nd - 00. (5) 

A powerful method for measuring 7 S in quasi-canonical simulations is provided by 
the size distribution of minbu's (minimal neck baby universes) jl2| — a part of the 
triangulation connected to the rest via a minimal neck. By counting in how many 
ways a minbu of size B can be connected to a surface of size (Nd — B), the size 
distribution becomes 

_ „ B Z B (N D - B)Z Nd _ b 

n No {B) k, (b) 

Zn d 

~ [(N D - B) B]^~ 2 . 

The distribution n^ D (B) is measured in simulations and 7 S determined by a fit to 
Eq. (6). 

(ii) The fractal or Hausdorff dimension dn measures the intrinsic "dimensionality" 
of the triangulations and is defined by the volume of a geodesic ball with radius r: 
v(r) ~ r du . The geodesic distance is defined as the shortest path dij, between two 
vertexes, traversed along links, either on the triangulations or its dual graph. 

To measure the fractal dimension dn one explores the scaling behavior of the 
vertex-vertex, or two-point, correlation function, 

1 



which counts the number of vertexes at distance r from a marked vertex i. Assuming 
that the only relevant length-scale in the model is defined by N D , general scaling 



arguments |13| imply that 

9N D (r) - N d - 1 ^f("). (8) 



N Vd H 



The optimal scaling of distributions gN D (r), corresponding to different volumes, 
onto a single scaling curve defines du- 
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1.3 The pure gravity phase structure 

Extensive numerical simulations have established that the model Eq. (|j) has two 
distinct phases, both in three and four dimensions. There is a strong-coupling 
(small k) crumpled phase and a weak-coupling (large k) elongated phase. The two 
phases are separated by a phase transition at critical value of the inverse Newton's 
constant k c . 

For k < k c the internal structure collapses and the crumpled phase is charac- 
terized by a singular structure, i.e. a small set of sub-simplexes connected to an 



extensive fraction of the total volume [ 14 ] . In D = 4 this structure consists of two 



singular vertexes, with local volumes lo' v that grows linearly with the volume N4, 

joined by a sub-singular edge; to' e ~ A^ 3 . In D = 3 there is only a sub-singular 
1 

vertex with oj' v ~ iV 3 . The geodesic distance between any two vertexes always 
stays at the level of the lattice cut-off. This prevents any sensible continuum limit 
in this phase — the quantum geometry is essentially zero-dimensional. In addition, 
large geometric out-growths are suppressed and, at least formally^], 7 S = — oo. 

In the elongated phase the internal geometry is composed of bubbles glued to- 
gether via small necks into tree-like structures — branched polymers — with du = 2 
and 7 S = 1/2. The neck size does not grow with the volume and the internal 
geometry is essentially one-dimensional. 

As the fractal structures described above are somewhat pathological, it is unlikely 
that an interesting continuum limit exist in either of the two phases. The initial hope 
was, however, that the transition point k c would correspond to a non-trivial fixed 
point where a sensible continuum limit could be defined. However, the transition 
turned out to be discontinuous, both in three ]l5| , [lG| ] and four |17]] dimensions; this 
is demonstrated in Figure 3. 



2 In the crumpled phase the sub-leading asymptotic behavior is Zn d (^) ~ exp(/iiVzj + aNp), 



b < 1, rather than Eq. Hence 7 S is not well defined. 
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Figure 3: Histogram of the energy density, no = Nq/Nd, for pure simplicial gravity at 
the transition k c . Both in (a) three |l6||, and (b) four jl7| dimensions a two state signal is 
observed, characteristic of a discontinuous phase transition. 



In four dimensions the transition is weakly first order and, consequently, hard 
to observe — it is only manifest on large volumes, N4 > 32.000. As present day 
computers limit practical simulations to N4 < 10 5 , this poses a serious problem for 
any numerical investigation of the model Eq. (|3|). 
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2 Modified Models of Simplical Gravity 



The discontinuous nature of the phase transition implies that the model Eq. (||) does 
not posses a continuum limit where a theory of gravity could be defined. Does this 
mean that there is something inherently wrong with the dynamical triangulation 
approach, or is the discretization Eq. (|3|) simply too naive — some fine-tuning of 
additional parameters might be needed. The success of two-dimensional models of 
dynamical triangulations, where the critical behavior is remarkable insensitive to 
details of the discretization, might be misleading. The theory of gravity in higher 
dimensions is more complicated and more carefully defined discretization may be 
necessary. 

This has led to investigations of modified models of simplicial gravity; modified 
either by changing the relative weight of the triangulations with e.g. a measure term, 
or by coupling matter fields to the theory. Although preliminary investigations of 
such modified models did not observe any qualitative change in the critical behavior 
[18, 19], more recent numerical simulations |2(], [TJ], 21] have revealed that such mod- 



ifications do indeed influence the phase structure: The phase transition appears to 
soften, or even disappear. Moreover, a new crinkled phase is observed, characterized 
by internal geometry closer to what one might expect of smooth manifolds. 

Two such modified models of simplicial gravity have recently been investigated 
in D > 2: 

(a) A modified measure: 

z Nd ( K ,p) = e n °("*)^ -OTp (9) 



where o(rii) is the local volume of an sub-simples i — the number of D-simplexes 
containing i. 



(b) Coupling to f copies of matter fields: 

i 

~C(T) 



Zn d (kJ) = £ ( Z m(T)) / (10) 



where Zm is the matter partition function. 

In Ref. [pifl the particular case of vector fields coupled to four-dimensional sim- 
plicial gravity was investigated: 

Zm(T) = /il ^(Oe-SvF (11) 

where A(l) are non-compact U(l) gauge fields residing on the edges I, A(l a b) = —A(lf, a ) 
and the sum is over all triangles t a b c in the four-dimensional triangulation. 
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Figure 4: The phase structure of two-dimensional simplicial gravity coupled to matter helds 
and/or a measure term. 



2.1 Lessons from two dimensions 



Before discussing the phase structure of the models Eqs. (Q) and (10) for D > 2, it 
is worth considering the impact the two modifications have on simplicial gravity in 
two dimensions.^] 

In two dimensions modifying the measure has no effect on the critical behavior, 
except for large negative values of the coupling (3 where formations of large curvature 
defects is enhanced and there is a cross-over to a crumpled phase, dominated by 
intrinsically collapsed configurations. 

In contrast, adding matter fields changes the critical behavior. For matter fields 
with central charge c < 1 this corresponds to coupling conformally invariant matter 
to 2D-gravity (Liouville theory). In this phase the fractal structure displays a 
remarkable degree of universality; it only depends on the total central charge of 
the matter fields: — oo < 7 s (c) < 0, and 2 < djy(c)<4. For c > 1, on the other 
hand, the fractal structure degenerates into branched polymers. This is understood 



qualitatively as a condensation of spiky configurations [22]; the free energy of spikes, 
-^spike ~ (1 — c ) l°g(l/ a )> estimated from the dynamics of the conformal factor, shows 
that such singular configurations dominate the partition function for c > 1. 

What does this imply for the phase diagram in higher dimensions? Is it possible 
that the similarity between the phase structure for pure gravity in D > 2 and the 
c > 1 region in Figure 4 is not coincidental. That a suitable tuning of the parameters 
(5 or /, could lead to some higher-dimensional analog of the Liouville phase? 

This scenario is supported by recent analysis in Ref. [23] which suggest that 
something analogous to the "c = 1 barrier" might exist in four dimensions. Using 



3 In D = 2 the phase diagram simplifies as the Einstein term in the action is a topological 
invariant and the dynamics is insensitive to the corresponding coupling constant (as long as the 
topology of the manifold is fixed). 
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Figure 5: A schematic phase diagram of modified models of simplicial gravity for D > 2. 



an effective action for the conformal factor, they calculate the free energy of spikes: 

Fspikc ~ ^(Q 2 -8) log(i). (12) 

Q 2 = m ( N s + T N F + 62N V ~ 28 ) + ^Grav P la y s tlie role of a central charge and 
Ns,f,v refers to the number of massless scalar, fermion and vector fields contributing 
to the trace anomaly. A one-loop calculation of the contribution from the transverse 
gravitons, Qq I£lv ~ 7.9, suggests that in the absence of matter fields the theory has 
no sensible stable vacuum. It is intriguing, however, that the contribution of matter 
and ghost fields appears in D = 4 opposite to that found in D = 2; in four dimensions 
adding matter fields might actually stabilize the theory! 



2.2 Modified phase structure in D > 2 

Motivated by these considerations, extensive numerical investigations have been 
carried out on the models Eqs. (||) and ( |T^ ) ]i~6| , pC| , 21]. In addition to MC simula- 
tions, a method for calculating the strong-coupling expansion of simplicial gravity in 
D > 2 has been developed pi] ], and has proved invaluable for a qualitative view of 
the impact different modifications have on the critical behavior. The results of these 
investigations, which suggest a phase diagram akin to Figure 5, can be summarized 
as follows: 

(i) As several vector fields are coupled to simplicial gravity (or the curvature 
fluctuation are enhanced by modifying the measure) the discontinuous phase 
transition (solid line), separating the crumpled and the elongated phases, is 
softened and eventually, at some value / > 3 (/3 < — 1), either becomes contin- 
uous or disappears altogether. 

(ii) At the same value of / (/3) the polymerizations of the geometry is sup- 
pressed, the elongated phase disappears, and a new crinkled phase appears. 
The crinkled phase is separated from the crumpled and elongated phases by 
lines of either continuous phase transitions or by a crossover region (dashed lines). 
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Figure 6: The exponent j s versus a modified measure (solid line) or coupling to matter 
fields (dashed) from the strong-coupling expansion for k 3> k c . The circles are from MC 
simulations of the model Eq. ph for k — 4.5 pi|. 



In contrast to two dimensions, a remarkable degree of universality appears in 
simplicial gravity in D > 2; all "reasonable" modifications seem to lead to the same 
qualitative phase structure Figure 5 [21]. 

The nature of the crinkled phase has been explored for D = 4, both in MC 
simulations and from the strong-coupling expansion plj . As the number of vector 
fields / is increased (or j3 decreased), the string susceptibility exponent j s decreases 
(Figure 6). Moreover, modulo some trivial re-scaling: (3 = — A/ — |, the results for 
the two models are identical. Measurements of the fractal dimension in the crinkled 
phase yield dn ~ 4. 

At the transition lines separating the crinkled phase from the rest of the phase 
diagram (the dashed lines), no divergence is observed in the specific heat, only a 
cusp. It is thus difficult to determine in MC simulations if these are soft continuous 
transitions (third or higher order) or if this represents a smooth cross-over to dif- 
ferent a fractal structure. However, a soft continuous phase transition would not be 
unexpected; coupling to geometric disorder tends to soften phases transitions. For 
example, most transitions for matter coupled to 2L>-gravity are of third order. 

How should one interpret the phase structure Figure 5? Does the crinkled phase 
play a role analogous to the Liouville phase in two dimensions, or is the internal ge- 
ometry just as pathological in this phase as in the crumpled and elongated phases? 
Alternatively, is a continuum limit possible at some of the observed transition lines; 
provided, of course, they correspond continuous phase transitions? Further investi- 
gations will hopefully shed some light on these issues. 
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2.3 The strong-coupling expansion 

The strong-coupling expansion — an explicit enumeration of all triangulations up 
to a given volume — provides a powerful tool for exploring the phase structure 
of different modifications of simplicial gravity. It has previously been used in two 
dimensions, where a matrix- model formulation allows a recursive construction of all 
triangulations [24]. Combined with an appropriate series extrapolation method, the 



strong coupling expansion yields, for some models, very accurate estimates of the 
critical exponents. 

In D > 2 a recursive algorithm for generating all triangulations is not available. 
Instead an alternative procedure constructing the series in numerical simulations 
has been developed |^l], |25|l : 

(i) Identify all distinct triangulations T using MC simulations and a suffi- 
ciently complicated "hash" function. 

(ii) Calculate the corresponding symmetry factor C(T) by comparing all per- 
mutations of the vertex labels. 

(iii) For each triangulation T calculate the appropriate contribution from 
matter fields and/or a modified measure. 

The hash function f{T) plays a vital role in the identification. It should be 
sufficiently complicated to distinguish between combinatorially inequivalent trian- 
gulations, yet simple enough to be repeatedly calculated. The identification is ver- 
ified by comparing the calculated symmetry factors C(T) to the relative frequency 
with which different triangulation appear in the MC simulations. In practice, the 
identification procedure is much simplified by observing that only triangulations not 
reducible by any of the volume decreasing geometric moves have to be considered 
in the MC; all others are systematically constructed from smaller ones. In practice 
99.9% of the triangulations turn out to be reducible. 

This method has been used to calculate the first 15 to 20 terms in the strong- 
coupling expansion of the partition function Eq. (|3|), both in three and four dimen- 
sions. This corresponds to roughly 10 6 distinct triangulations. Given the triangu- 
lations it is easy to calculate the contribution from different types of matter fields, 
or from a modified measure, to the partition function. To determine the critical 
behavior, appropriate series extrapolation methods are used, e.g. the ratio method 



[24], in which higher order corrections to the partition function are systematically 
eliminated. 

The method is, however, limited to the region of the phase-space where the par- 
tition function is rapidly convergent. In cases where the sub- leading corrections to 
the partition function are not analytic, such as in the presence of a confluent singu- 
larity, the extrapolation fails to converge. In simplicial gravity the strong-coupling 
expansion has provided valuable informations about the nature of the elongated and 
crinkled phases; however, it has been of limited use in the crumpled phase and in 
the transitions region. 
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Figure 7: The phase diagram of balls-in-boxes 
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2.4 "Balls in boxes" 

Many of the characteristics of the phase structure of simplicial gravity, described 
in the previous sections, are captured by a simple mean-field model — the balls-in- 
boxes model [S3] . This model consists of fixed number N of balls distributed into a 
variable number M of boxes, 1 < M < M max . In its partition function, 



z n {k,p) = eKm E K<?l)-KSm) s, 



qi + ...+q m ,N i 



(13) 



m=l 



(li- 



the partitions of balls are weighted by the product of one-box weights p(qi) = q~^ , 
where q\ = 1, 2, ... is the number of balls in box i. This particular choice of weights 
mimics the modified measure, Eq. (||), introduced in simplicial gravity. 

The model Eq. (13) derives its non-trivial features from the constraint on the 
total number of balls. Introducing the average density of balls p = M/N, the "cur- 
vature", there is phase transition at a critical density p c . This transition separates 
an elongated (fluid) phase, with the balls evenly distribution among the boxes, from 
a crumpled phase, where the balls "condensate" into one singular box (akin to Bose- 
Einstein condensation). The transition is first order for (3 £ (2, oo), whereas it is 
continuous for /3 6 (1,2]. 

By introducing an additional kinematic bound on the average density of balls 
per box, p < p art , a new condensed phase appears in the model for (3 > 2.5 [27]. This 



is shown in Figure 7. In this phase, which resembles the crinkled phase of simplical 
gravity, the model is neither crumpled nor fluid. The average density per box is 
finite, (r) > 0; however, there is a singular box in the system which captures a finite 
number of the balls: {q S mg)/N > 0. 
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It is remarkable how such a simple model successfully describes many of the 
geometric features of simplicial gravity in D > 2, e.g. the appearance of a singu- 
lar vertex and the phase structure, while ignoring all the non-trivial information 
encoded in the connectivity of the triangulations. 



3 Degenerate Triangulations 

In the partition function Eq. (0), T denotes a suitable ensemble of triangulations 
included in the sum. Different ensembles are defined by imposing various restriction 
on how the simplexes are glued together. Provided this leads to a well-defined parti- 
tion function, and as long as this difference is only at the level of the discretization, 
one expects different choices of T to lead to the same continuum theory in the ther- 
modynamic limit. This is, of course, not true if the partition function is divergent 
as is the case even in two dimensions if the topology is not restricted. 

Traditionally, simulations of simplicial gravity in D > 2 have used an ensemble of 
combinatorial triangulations Tq- In a combinatorial triangulation every D-simplex 
is uniquely defined by a set of D + 1 distinct vertexes — it is combinatorially unique. 
This implies that two distinct simplexes have at most one face in common; this 
amounts to exclude both tadpoles and self-energy diagrams in the corresponding 
dual graph. Easing this restriction, allowing either multiple connected simplexes 
or simplexes glued onto themselves, defines ensembles of restricted and maximally 
degenerate triangulations, Td r and 1b M , respectively |3(J]. In the former ensemble 
to each simplex there is associated a set of D + 1 distinct vertexes; however, the same 
set of vertexes may be shared by any number of distinct simplexes. In contrast, in a 
maximally degenerate triangulation even the set of vertexes associated to a simplex 
may be degenerate, i.e. it may contain the same vertex more than once. Clearly, 
Tq C Td r C Td m ■ Figure 8 shows two-dimensional examples of those three different 
types of triangulations. 

In two dimensions the model Eq. (|3|) defined with those different ensembles is 
solvable as matrix model and, in all three cases, it yields the same continuum theory 
[28]. Then why consider different ensembles; why not use the one most convenient 
in the simulations? It turns out that the finite-size effects depend strongly on the 



ensemble used [29|. In particular, the finite-size effects are considerable reduced 
the less restrictions are placed on the triangulations. An example of this is the 
volume dependence of the effective magnetization exponent, P/vdn, for an Ising 
model coupled to 2D-gravity (Figure 9). This reduction in the finite-size effects is 
understandable — with a large ensemble of triangulations it is easier to approximate 
a given continuum geometry using triangulations of finite volume. 



3.1 D>2 

It is only recently that simulations of simplicial gravity in higher dimensions have 
been extended to include degenerate triangulations. In three dimensions both re- 
stricted and maximally degenerate triangulations have been investigated fl3Cf| ; how- 
ever, only the former ensemble was found to lead to a well-defined partition function 
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Figure 8: Three types of two-dimensional triangulations: (a) combinatorial, (6) restricted 
degenerate, and (c) maximally degenerate. The dashed lines indicate the corresponding 
dual graphs. 
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Figure 9: The effective magnetization exponent, /3/i/du, for an Ising model coupled to 
2D-gravity versus volume p9[ |, both for combinatorial (squares) and degenerate (circles) 
triangulations |29 . The theoretic value is 1/6. 
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Figure 10: The pseudo-critical cosmological constant fi c (N4) for ensembles of AD com- 
binatorial (squares) and degenerate (circles) triangulations. Included are fits assuming an 
power-law convergence (solid lines) and a logarithmic divergence (dashed lines) |3ll . 
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Eq. (|3|). More recently restricted degenerate triangulations have also been used in 
four dimensions 

Qualitatively, the phase structure does not change when degenerate triangula- 
tions are included in the model Eq. (|3|). There is still a small k crumpled phase 
and a large k elongated phase, separated by a discontinuous transition. However, 
just as in two dimensions the finite-size effects appear much reduced, especially in 
the crumpled phase. An example of this is the pseudo-critical cosmological constant 
H c (Nd), which converges much more rapidly for the ensemble To fl than for Tq- This 
is shown in Figure 10 for D = 4 and k = 0. 

This important observation demonstrates (numerically) an exponential bound on 
the number of degenerate 4D triangulations of fixed volume^ - a necessary condi- 
tion for the convergence of the partition function Eq. (||). This is demonstrated in 
Figure 10 by two fits: one fit assumes an exponential bound on the canonical par- 
tition function (solid), the other a super-exponential growth (dashed). Whereas for 
combinatorial triangulations the fits are somewhat inconclusive, for degenerate tri- 
angulations a super-exponential growth is definitely ruled out [31|. And, as Tq C Tjy, 
this automatically implies a bound on the entropy of combinatorial triangulations 
as well. 



4 Physical Membranes 

The study of the elastic and geometrical properties of physical membranes, two- 
dimensional surfaces with an extrinsic curvature embedded in R 3 , is a rapidly ad- 
vancing field with a nice interplay between experiments, analytical theory and com- 



puter simulations (see e.g. Refs. [33, |34j). Viewed as two-dimensional generalization 
of one-dimensional chains, membranes show a rich behavior on mesoscopic length 
scales. They can be classified according to their internal structure: 

(a) Polymerized or crystalline membranes have fixed internal structure and are 
characterized by both long-range orientational and translational order. 

(b) Fluid membranes, with dynamical internal connectivity and no internal 
order. 

An example of the former is the spectrin cytoskeleton of red bloods cells, whereas 
bilayers of amphiphiles such as phospholipids generally exhibit liquid-like behavior. 
Physical membranes encountered in nature are, of course, self-avoiding; a property 
pivotal to many aspects of their phase structure. They are, however, commonly 
studied as self-intersecting (phantom) membranes, both for simplicity and as statis- 
tical system interesting in their own right. 

It is by now well established, both by analytical and computational work, that 
crystalline membranes exhibit a low-temperature flat (ordered) phase and, at least 
for phantom membranes, a high-temperature crumpled phase [j35|. The existence 
of an ordered phase in a two-dimensional system with a continuous symmetry and 
an (apparent) local interaction is remarkable given the Mermin-Wagner theorem. 

4 The existence of an exponential bound has been somewhat controversial for the ensemble of 
combinatorial triangulations in four dimensions [BJ. 
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What stabilizes the flat phase in a crystalline membrane are the out-of-plane fluc- 
tuations that couple to the in-plane "phonon" degrees of freedom — bending of the 
membrane is accompanied by internal stretching. 

The same argument for a stable flat phase does not apply in the case of a 
fluid membrane where the internal stretching can be compensated by a flow of 
"particles" into the distorted area, thereby "screening" the curvature fluctuations 
from the rest of the membrane. The absence of a stable flat phase is supported by 
renormalization group analysis of models of fluid membranes, which suggest that 
at large distances the bending rigidity becomes irrelevant [pq| . This is, however, 
contradicted by numerous numerical simulations of fluid membranes that indicate a 
crumpling transition, akin to that of its crystalline counterpart [p^] . 



4.1 Including anisotropy 

The above discussion applies to isotropic membranes. Remarkable, adding intrinsic 
anisotropy to a model of crystalline membrane profoundly influences the global phase 



diagram [38|. Most notably, a new tubular phase appears, characterized by the 
presence of long-range orientational order in one direction only — in the transverse 
directions the membrane is crumpled. For a phantom membrane, the tubular phase 
separates the existing crumpled and flat phases; only an isotropic membrane passes 
directly from a flat to a crumpled state (Figure 11). Analogous to a flat membrane, 
a tubular membrane is stabilized by the transverse stretching energy cost of bending 
fluctuations in the extended direction. 

The existence of a tubular phase has been confirmed in recent numerical simula- 
tions^] on a crystalline membranes |39| . The membrane is discretized on a triangular 

5 Simulations of physical membranes are notoriously difficult due to very long auto-correlation 
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mesh, with the topology of a disk, and the interactions modeled by Hamiltonian 
composed of two terms: a Gaussian pair potential between neighboring nodes, and 
a bending energy introduced as ferromagnetic interaction between normals n a on 
adjacent triangles: 

= ^2 l r ^~ r ^'| 2 ~ • n 6 - K 2 ^ (y) n a • n 6 . (14) 

(era') (ab) (ab) 

The anisotropy is introduced as different bending rigidity in the two intrinsic direc- 
tions, x and y. 

The phase diagram of the model Eq. (|14|) was explored, for different bending 
rigidities (ki , K2), and evidence for two distinct phase transition — a crumpled-to- 
tubular and tubular-to-flat — is observed in the fluctuations in the two bending 
energy terms: C'y{n) and C v {n) (Figure 12). 

In the tubular phase, the tubule cross-section radius R± (the radius of gyration) 
and the undulations /i rms transverse to its average axis of orientation, scale like: 
~ L Up and h rms ~ L^, where vp and £ are the size (Flory) and roughness 
exponents. MC simulations of a tubular membrane give: vp = 0.269(7) and £ = 
0.850(40) |39|], compared to the theoretical predictions: up = -? and £ = 1 [38]. 



Notice that £ = 1 implies that a phantom tubule is only marginally stable with 
respect to the fluctuations in the transverse directions; the tubule almost crumples. 



times present in standard MC simulations. These auto-correlations are, though, much reduced 
by applying more powerful updating methods such as overrelaxation or unigrid algorithms in the 
simulations |40t|. 
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4.2 Self-avoiding membranes 



The more physical self- avoiding membranes present more of a challenge, especially 
for numerical simulations, due to the non-local nature of the self-avoiding constraint. 
For a flat isotropic crystalline membrane self-avoidance is expected to be irrelevant; 
self-avoidance is, however, expected to suppress the crumpled phase. Hence only 
flat self-avoiding isotropic membranes are expected to exist; this is supported by 
numerical simulations. 

For membranes crumpled in one direction only, self-avoidance is less constraining 
and the tubular phase is expected to survive in the more physical self-avoiding case 
[38, 41, 42]. However, unlike for a flat membrane self-avoidance is expected to be 
relevant for tubules. Adding a non- vanishing self- avoiding coupling leads to a flow 
in the couplings to a new critical (infrared stable) fixed point, distinct from the one 
that governs the critical behavior of a phantom tubule. This fixed point structure 



is supported by renormalization group analysis [|38| , 42]. An improved one-loop 
calculations of the critical exponents v and in an e-expansion around the upper 
critical dimension dc P = 11, gives: v = 0.62 and £ = 0.8 [42]. 

It is of great interest, albeit very challenging, to verify this critical behavior in 
numerical simulations of a self-avoiding anisotropic membrane. This work is under- 
way [43]. Preliminary results indicate the existence of only one phase transition, 
from a tubular to a flat phase; as expected a crumpled phase is not observed. 



5 Outlook 

(i) Much work is needed before the phase structure of modified models of simplicial 
gravity in D > 2 is fully understood. The nature of the crinkled phase, and of the 
observed phase transitions, must be established before any statements can be made 
about the relevance the observed phase structure Figure 5 has for any potential 
continuum theory of quantum gravity. 

(ii) Degenerate triangulations might prove a valuable tool in any such investiga- 
tions. In addition to reducing the finite-size effects in simulations, they also serve 
as a consistency check on any observed phase structure in simplicial gravity — any 
critical behavior that is relevant in a continuum limit should be independent of the 
ensemble of triangulations used in defining the model. 

(iii) Anisotropic crystalline membrane have a surprisingly rich phase structure, 
some of which has been revealed through analytic calculations and/or numerical 
simulations. To extend these simulations to the more physical self-avoiding case is 
of much interest, albeit very challenging. 

Fluid membranes present another challenge to numerical simulations. They 
are expected always to crumple on long length-scales, yet computer simulations 
suggest a possible transition to a flat phase. However, these results where obtained 
is simulations of systems of modest size; revisiting fluid membranes now that better 
algorithms and faster computers are available would be an interesting project. 
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